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Показано что применение ортогонализации Г оама-Шмидта векторов измеренных значений при рекуррентном оценивании па- 
раметров модели объекта управления позволяет сократить количество арифметических операций за счет отказа от процедуры 
псевдообращения получаемой матрицы. На каждой итерации алгоритма оценивания необходима ортогонализация только од- 
ного текущего вектора измеренных значений. Такой подход приводит к существенному повышению быстродействия алгоритма 
оценивания. 


1. Введение 

Как известно [1,2], наряду со скоростью сходи- 
мости и помехоустойчивостью важнейшим свой- 
ством рекуррентных алгоритмов оценивания пара- 
метров линейных и нелинейных моделей объектов 
управления является их экономичность, характе- 
ризующаяся количеством арифметических опера- 
ций над вещественными числами, необходимых 
для реализации одной итерации уточнения оценок 
параметров идентифицируемых моделей. Особен- 
но актуальным требование экономичности данных 
алгоритмов оказывается в тех случаях, когда выде- 
ляемое на их реализацию время в одном такте 
функционирования системы управления является 
в значительной мере ограниченным. С подобного 
рода ситуациями неизбежно приходится сталки- 
ваться, например, в тех случаях, когда управляе- 
мым объектом является некоторый быстропроте- 
кающий процесс, для эффективного управления 
которым требуется достаточно высокая частота 
контроля его состояний и коррекция управляющих 
воздействий. Аналогичные ситуации возникают и в 
тех случаях, когда управляемый объект не является 
каким-либо быстропротекающим процессом, но 
для коррекции управляющих воздействий требу- 
ются значительные затраты машинного времени. 
Последнее же, очевидно, может иметь место из-за 
недостаточно высокого быстродействия использу- 
емой в системе ЭВМ, либо из-за значительных 
объемов вычислений, необходимых для нахожде- 
ния новых управляющих воздействий и других за- 
дач, решаемых системой управления. Обе эти си- 


туации являются достаточно типичными при раз- 
работке автоматизированных систем управления, 
базирующихся на мини- и микроЭВМ, быстродей- 
ствие которых, как правило, не столь велико, как 
это часто бывает необходимо. 

В данной работе предлагается модификация 
многоточечного рекуррентного алгоритма оцени- 
вания параметров моделей линейных статических 
объектов управления, основанного на использова- 
нии псевдообратных матриц, требующая для своей 
реализации существенно меньших объемов вычи- 
слений, чем это необходимо для реализации моди- 
фицируемого алгоритма. Значительное сокраще- 
ние объемов вычислений здесь удается добиться за 
счет ортогонализации измерений входных пере- 
менных объекта с применением хорошо известной 
процедуры Грама- Шмидта ортогонализации векто- 
ров [3, 4]. 


2. Постановка задачи рекуррентного оценивания 
параметров математических моделей объекта 
и описание алгоритма ее решения, основанного 
на использовании псевдообратных матриц 

Задачу рекуррентного оценивания параметров 
модели управляемого объекта сформулируем сле- 
дующим образом. Пусть имеется линейный стати- 
ческий объект, значения «-мерного входа и скаляр- 
ного выхода которого связаны соотношением 


У, = 


4 . \ 

х,,а 

у 


( 1 ) 
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где х ( = (х ІІ ,х і2 ,...,х а ), у - соответственно значе- 
ния входа х и выхода у в момент времени /; 

і т 

а = (а 1 ,а 2 ,...,а п ) - л-мерный вектор неизвестных 
параметров; Т - символ транспонирования векто- 
ров и матриц; я - конечное натуральное число. 
Пусть, далее, в каждый момент времени / имеются 
измерения вида 



л 

Х г 



Х,= 

Я-1 

4 

и У, = 

У ,- 1 


^,-му 


\Уі-і + 1 у 


В настоящее время существует целый ряд алго- 
ритмов, позволяющих осуществлять ее решение 
[2]. Одним из них является многоточечный рекур- 
рентный алгоритм оценивания параметров линей- 
ных моделей (1), основанный на применении псев- 
дообратных матриц [5], который имеет вид 

'I' >1' + >1' 4'' 

а ,=а 1 _ ] + Х і (у І -Х і я,_, ), (3) 

где X, - (/хя)-матрица, псевдообратная к матрице 
X,. Алгоритм вычисления данной матрицы в каж- 
дый дискретный момент времени / составляют сле- 
дующие операции: 

X, =х* = х 1 (х 1 ,х і у', (4) 


4 

где Х г у ( - соответственно (/хя)-матрица и /-мер- 
ный вектор измерений входа х и выхода у объекта; 
/ — глубина памяти алгоритма - некоторое нату- 
ральное число, такое, что 1 </<Л/; М— конечное на- 
туральное число, которое может быть как больше, 
так и меньше числа л. Задача рекуррентного оцени- 
вания параметров модели объекта при каждом / за- 

4. 

ключается в том, чтобы получить новые оценки я, 

параметров а модели ( 1 ) на основе измерений 

4- 

(2) и уже имеющихся к моменту / оценок я,_, дан- 
ных параметров, полученных на предшествующем 
(/ — 1 )-м такте оценивания. 

Приведенная постановка задачи рекуррентного 
оценивания параметров модели является одной из 
простейших задач идентификации математических 
моделей управляемых объектов. Вместе с тем с точ- 
ки зрения практических приложений именно она 
имеет важнейшее значение и представляет перво- 
степенный интерес. Это обуславливается, во-пер- 
вых, тем, что количественные связи между пере- 
менными многих реальных объектов имеют линей- 
ный или весьма близкий к нему характер и, следо- 
вательно, с достаточно высокой точностью могут 
быть описаны моделями вида (1). Во-вторых, доста- 
точно типичной является ситуация, когда исследуе- 
мым объектом является реальный технологический 
процесс, значения переменных которого изменя- 
ются в достаточно узких пределах, и поэтому, если 
под значениями у, и х, понимать не значения пере- 
менных у их, а их отклонения от некоторых фикси- 
рованных, например, номинальных или средних 
значений, то количественные связи между данны- 
ми отклонениями можно, очевидно, также с доста- 
точно высокой точностью описать моделями вида 
( 1 ), несмотря на заведомо нелинейные зависимости 
между значениями переменных у и х. И, наконец, 
в-третьих, с помощью тех или иных взаимноодноз- 
начных функциональных преобразований пере- 
менных у их достаточно часто удается перейти к 
переменным у' и х\ количественные связи между 
которыми адекватно описываются моделями (1). 


Х к = 


X, 


(х к -і к х *-і) + (& ф \ х к-і 

(1 + (? 4 ,г,))" 1 А' і _ 1 г/(г, |-1)Д =\ Х к 


( 5 ) 

к = \,1 


где х к - к- я строка матрицы X,; Х к - матрица раз- 
мерности (кхц), составленная из первых к строк 
матрицы X;, Х к - псевдообратная к Х к матрица, 

_> _> + ^ 

г = х к Х к _ } , 0„ - нулевой я-мерныи вектор. 


Основные свойства алгоритма (3) изложены в 
[2, 6] и сводятся к следующему. 

1) Глубина памяти I может быть меньше или равна 

4. 

размерности я вектора параметров а модели ( 1 ) 
либо больше ее. Данное свойство является след- 
ствием того, что псевдообратная матрица X, су- 
ществует при произвольных соотношениях 
между / и я, что открывает широкие возможно- 
сти выбора глубины памяти /. Выбирая глубину 
памяти, можно обеспечить эффективное при- 
менение алгоритма. 

2) Необходимым и достаточным условием сходи- 

4 

мости оценок я, , /=1,2,3..., вычисляемых с по- 
мощью алгоритма (3), является линейная неза- 
висимость измерения х, от предшествующих 
1-1 измерений входа и выхода объекта. 

4 

3) Последовательность оценок я, , вычисляемых с 
помощью данного алгоритма, является моно- 
тонно по евклидовой норме сходящейся к ис- 

4 

тинным значениям параметров а . 

4 

4) Сходимость оценок я, имеет место при произ- 

4 4 

вольных начальных оценках я 0 параметров а . 


5) Скорость сходимости алгоритма (3) монотонно 
увеличивается с увеличением глубины памяти и 
ее изменении в пределах от /= 1 до /=л. 
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6) Имеются возможности повышения помехоу- 
стойчивости рассматриваемого алгоритма за 
счет регуляризации алгоритма (4), (5) вычисле- 
ния матриц X, [2,7]. 

7) Количество арифметических операций, выпол- 

і 

няемых на одном такте уточнения оценок а п 
определяется следующими соотношениями: 

а)2Ѵ е =0,5(Зи-1Х/+1), 

б) N у = 0,5и(3/ + 5)1 и в) Ы 6 = I, (6) 

где /Ѵ с , Ы у и ЛГ,- соответственно, количество опера- 
ций сложения, умножения и деления веществен- 
ных чисел. 

3. Модификация многоточечного рекуррентного 
алгоритма (3) на основе использования 
ортогонализации измерений 

Как показано в работах [2, 6], трудоемкость ре- 

4 - 

ализации одного такта подстройки оценок а,_ х па- 

і 

раметров а модели (1) в соответствии с многоточеч- 
ным алгоритмом (3) характеризуется соотношения- 
ми (6). Непосредственный подсчет количества ариф- 
метических операций, необходимых для реализации 
формул (4)— (5) позволяет получить соотношения 
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! 

' х, л 
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и воспользоваться правилами псевдообращения 
блочных матриц. 

Вторая из обсуждаемых возможностей, реали- 
зованная в данной работе, состоит в том, чтобы 
воспользоваться каким-либо способом ортогона- 
лизации строк матрицы Х„ 1= 1,2,3,... Для этих це- 
лей, в частности, можно воспользоваться извест- 
ной процедурой Грама-Шмидта ортогонализации 
заданной совокупности векторов. В результате 
применения данной процедуры к строкам матрицы 
X, и объединения ее с алгоритмом (3) получаем мо- 
дифицированный многоточечный алгоритм, 

имеющий следующий вид: 

Г _ і _ _ _ 

ѵ, = х, + ^(х, , н>і )и’,-і, (8) 

1=1 

УѴі = ѴіІ\\ѵі II, (9) 

Ѳ,=[у,~Іі(х„т- 0^_ 1 ] / /||ѵ,||, (10) 

а, =а,_ 1 + тЩ-{ и-, ,а,- Х )\, / = 1,2,..., (11) 


а)7Ѵ с = 0, 5((3/ — 1)/і — / — 1)7, 
б) Ы у = 0, 5и(3/ + 1)/ и в) Ы д = /, (7) 

+ + + 

где символами И с , И ѵ , И г) обозначены соответ- 
ственно количество операций сложения, умноже- 
ния и деления вещественных чисел. 

Сопоставляя полученные соотношения с (6), 
можно видеть, что подавляющая часть вычислений 

г 

при вычислении оценок я, в соответствии с алго- 
ритмом (3) затрачивается на вычисление псевдооб- 
ратной матрицы X,- Вместе с тем, сравнивая измере- 
ния (2) на (/ — 1)-м и 1-м тактах оценивания, можно 
видеть, что только первая строка х, (компонента у,) 
г 

матрицы X, (вектора у,) является новой, а осталь- 
ные ее строки (его компоненты) являются строками 

і 

(компонентами) матрицы Х,_ х (вектора у м ). Отме- 
ченные обстоятельства наводят на мысль организо- 
вать вычисление X, так, чтобы в максимально воз- 
можной мере использовать результаты, полученные 
при вычислении Х,_ х , и за счет этого сократить 
объем вычислений на каждом такте оценивания па- 
раметров модели (1). Имеется по крайней мере два 
способа реализации данной идеи. Первый из них 
заключается в том, чтобы выделить в матрицах X, и 
Х,_ { блок строк, входящих одновременно в обе ма- 
трицы, т. е. представить их в виде 


где верхний предел суммирования Ь определяется 
соотношением 

Гг-1, 1 = 1,2,...,/— 1, 

[/- 1, 1 = 1 , 1 + 1 ,... . 

Здесь ||ѵ,|| - евклидова норма вектора ѵ„ вычисля- 


емая в соответствии с равенством ѵ< 

Из соотношений (8)-( 11) видно, что на первых 

г 

1-1 тактах оценивания параметров а осуществля- 
ется накопление измерений (2) и их ортонормиро- 
вание, а на каждом из последующих тактов - сдвиг 
/ ранее поступивших измерений с вычеркиванием 
самого «старого», т. е. (?— 1)-то измерения и орто- 
нормирование самого «нового» ?-го измерения по 
отношению к остальным 1-1+1 имеющимся изме- 
рениям. Поскольку сдвиг измерений и вычеркива- 
ние одного из них из ортонормированной совокуп- 
ности не нарушает ортонормированности остаю- 
щихся в ней измерений, то в обоих случаях необхо- 
димо ортонормировать только одно вновь посту- 
пающее измерение х,. Непосредственная коррек- 
і 

ция оценок я (-1 в этом случае также существенно 
упрощается. Общий объем вычислений, выполняе- 
мых на одном такте в установившемся режиме оце- 
нивания (при />/), в данном случае характеризуется 
следующими соотношениями: 




ѴН ) 
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Ы с = (21 + 1)и, /Ѵ ѵ = (3и + 1)/, /V., = 1 . (12) 

Сопоставляя данные соотношения с (6) и (7), мож- 
но видеть, что применение ортонормирования измере- 
ний (2) существенно сокращает объем вычислений. 
Так, зависимость величин 7Ѵ С и И у от глубины памяти / в 
(4) является квадратичной, а в ( 12) она имеет линейный 
характер. Значение же 7Ѵ Й в (12) вообще не зависит от 
значения /, а в соотношениях (6) значение /Ѵ„ равно /. 

Заметим, что сравнение экономичности алго- 
ритмов здесь производится на основе сопоставле- 
ния количества арифметических операций, выпол- 
нение которых необходимо для реализации одного 
такта подстройки параметров идентифицируемой 
модели в соответствии со сравниваемыми алгорит- 
мами. Подход к сравнению экономичности алго- 
ритмов на основе отмеченных выше критериев 
корректен лишь при использовании данных алго- 
ритмов в следящих режимах оценивания [1,2], ког- 
да скорости их сходимости примерно равны. Срав- 
нение же алгоритмов по данным критериям при 
использовании режима обучения не вполне кор- 
ректно, так как в этих режимах они имеют суще- 
ственно различные скорости сходимости. В подоб- 
ных ситуациях более корректным было бы сравне- 
ние экономичности алгоритмов по суммарному 
количеству арифметических операций, необходи- 
мых для реализации процесса подстройки параме- 
тров в целом, продолжающегося с первого такта 
оценивания и, например, до тех пор, пока не будет 
достигнута заранее заданная точность оценок. 

Однако использование такого подхода к оценке 
экономичности рекуррентных алгоритмов оцени- 
вания возможно лишь в некоторых частных слу- 
чаях и по результатам экспериментальных исследо- 
ваний в одних и тех же условиях. В общем же слу- 
чае его применение связано со значительными 
трудностями, так как скорость сходимости данных 
алгоритмов определяется не только их свойствами, 
но и свойствами измерений на каждом из тактов 
оценивания, и, прежде всего, теснотой линейной 
зависимости между ними и их точностью. 

Как видно из представленных выше результа- 
тов, использование ортогонализации измерений 
переменных объекта в значительной мере изменяет 
структуру алгоритма (3) и позволяет существенно 
повысить его экономичность. Изменяются при 
этом и некоторые другие свойства данного алгорит- 
ма. В частности, глубина памяти / алгоритма 
(8) — ( 1 1) не может быть больше п. Однако многие 
свойства алгоритма (3) при этом сохраняются и 
полностью наследуются алгоритмом (8)— (11). В 
частности, полностью сохраняются такие его важ- 
нейшие свойства, как сходимость при произволь- 

г 

ных начальных оценках я 0 и любом значении /, 
удовлетворяющем отмеченным выше ограниче- 
ниям, а также монотонность сходимости по евкли- 

4, і 

довой норме оценок я, к а при 1— »со. Справедли- 
вость всех этих свойств легко устанавливается с по- 


мощью тех же самых рассуждений и приемов, кото- 
рые были использованы при исследовании алгорит- 
ма (8) — ( 11) с увеличением его глубины памяти / в 
пределах от 1 до п. Справедливость данного свой- 
ства становится совершенно очевидной, если иметь 
в виду, что в случае ортогональных и точных изме- 
рений переменных объекта для получения вектора 
г 

а нужно использовать не более чем п измерений. 

Заключение 

Результаты многочисленных эксперименталь- 
ных исследований алгоритма (8) — (1 1), выполнен- 
ных в тех же условиях, что и исследования алгорит- 
ма (3), достаточно подробно рассмотрены в [7]. От- 
метим следующие основные выводы, вытекающие 
из данных результатов. Во-первых, они иллюстри- 
руют работоспособность исследуемого алгоритма и 
подтверждают наличие у него всех отмеченных вы- 
ше свойств. Во-вторых, при малых глубинах памя- 
ти (/< 5) алгоритм является устойчивым по отноше- 
нию к ошибкам измерений переменных объекта. 
В-третьих, при более значительных глубинах памя- 
ти (1> 6) на некоторых тактах оценивания наблюда- 
ются случаи неустойчивого поведения алгоритма, 
г 

когда оценки я, оказываются более удаленными от 

г і 

а , чем оценки я, . В-четвертых, неустойчивость 
алгоритма легко устраняется с помощью его регу- 
ляризации, осуществляемой заменой множителя 
||ѵ,||- 1 в формулах (9) и (10) множителем (||ѵ,||+г) _1 , где 
г- параметр регуляризации, что в вычислительном 
отношении незначительно усложняет весь алго- 
ритм. Наконец, в-пятых, наличие у алгоритма регу- 
лируемой глубины памяти / и параметра регуляри- 
зации г позволяет в каждом конкретном случае 
подбирать их значения таким образом, чтобы обес- 
печить желаемое быстродействие и помехоустой- 
чивость алгоритма. Проблема выбора оптимальных 
значений параметра регуляризации и глубины па- 
мяти алгоритма требует проведения дополнитель- 
ных теоретических и экспериментальных исследо- 
ваний и выходит за рамки настоящей статьи. Здесь 
мы отметим лишь следующее. Как можно меньшее 
значение глубины памяти / следует выбирать в слу- 
чае, когда необходимо обеспечить максимальное 
быстродействии и более высокую помехоустойчи- 
вость алгоритма. Значение же параметра регуляри- 
зации г на начальных тактах оценивания параме- 
тров модели управляемого объекта следует устана- 
вливать равным нулю или достаточно близким к 
нему. А на последующих этапах подстройки оценок 

і г 

я ( , т. е. при 1>п, когда оценки я, оказываются 

і 

близкими к истинным значениям параметров а, 
значение г должно быть отличным от нуля и тем в 
большей мере, чем большую помехоустойчивость 
необходимо обеспечить. 
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